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We have reconstructed the radiative decays Xbi^P) — > T(IS') 7 and 
XbC^P) — > T{1S) 7 in pp coUisions at ^/s = 1.8 TeV, and measured 
the fraction of T(IS') mesons that originate from these decays. For 
T(IS') mesons with pj. > 8.0 GeV/c, the fractions that come from 
X6(1P) and XbC^P) decays are {27.1±6.9{stat)±4A{sys))% and (10.5± 
4.4(stat) lb 1.4(sys))%, respectively. We have derived the fraction of 
directly produced T(IS') mesons to be (50.9 it 8.2{stat) ± 9.0{sys))%. 

PACS numbers: 13.85.Ni, 14.40. Gx 

The large discrepancies between the charnionium production cross sections mea- 
sured by the Collider Detector at Fermilab (CDF) |jl[ and the predictions of the 
Color Singlet Model (CSM) can be explained in a theoretical framework based on 
non-relativistic QCD 0]. In this model, originally developed to describe rigorously 
the decay of heavy quarkonium states, the production process is factorized into short 
distance cross sections to produce the heavy quark pair, and long distance matrix 
elements, describing their binding into the quarkonium state. These matrix elements 
must be determined from experimental data but are assumed to be independent of 
the reaction and can be used to predict other processes. A consequence of this ap- 
proach, when applied to charmonium production in pp collisions, is the realization 
that cc pairs, produced at short distance in a color-octet state, are responsible for 
the bulk of the cross section. In the bottomonium sector CDF has measured the in- 
clusive production cross section of T(IS'), T{2S) and T{3S). The prediction of CSM 
underestimates the measured rate, although by a smaller amount than found for char- 
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monium 0]. Color-octet contributions can account for the discrepancies, but data on 
the inclusive T cross section alone is not enough to extract the matrix elements with- 
out assumptions [Q. In order to do this one needs to separate experimentally T's 
produced directly from those arising from the decays of heavier mesons. 

In this letter we report a study of the reaction pp -^ XbX, Xb -^ '^{^S)-y, T(IS') -^ 
fi~^fi~ at ^/s = 1.8 TeV using CDF. This analysis, based on approximately 90 pb~^ of 
data collected during the 1994-1995 collider run, describes the first observation of xt 
mesons at a hadron collider. Since the branching fractions for xt decays into other 
modes containing an T(IS') are expected to be small, this study allows us to measure 
the contribution of Xb decays to T(15') production. Even though T mesons can be 
reconstructed at CDF throughout the low p^ region, we perform this measurement 
only for p"^ > 8.0 GeV/c because at lower p^ the photon emitted in the radiative Xb 
decay is not energetic enough to be detected efficiently. In this analysis we do not 
study transitions of Xb niesons to T{2S) because photons from this decay have even 
lower energy. 

The CDF detector has been described in detail elsewhere 0. The events used 
in this analysis were collected with a three-level trigger system which selects events 
consistent with the presence of two muons. The first level required that two candidates 
be observed in the muon chambers. The second level required that two or more 
charged particle tracks, partially reconstructed in the central tracking chamber (CTC) 
using a fast processor, matched within 15° in (p (the azimuthal angle) the muon 
candidates. The third level required better precision on the azimuthal matching and 
required the dimuon invariant mass to be between 8.5 and 11.4 GeV/c^. 
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To identify T's we select pairs of oppositely charged muon candidates with 
Pt > 2.0 GeV/c, and we require the pair to have pT(/U^/i^) > 8.0 GeV/c. Since 
T mesons do not originate from long lived particles we constrain the muon tracks 
to originate from the primary interaction point to improve mass resolution. Fig. [^ is 
the resulting dimuon invariant mass distribution showing three peaks corresponding 
to the T(15'), T{2S) and T{3S) resonances. Due to the trigger and muon acceptance, 
the pseudorapidity of the muon pairs is limited to the central region, corresponding ap- 
proximately to |r7(/i+/i^)| < 0.7, where rj = — ln[tan(^/2)] and 6 is the polar angle with 
respect to the beam axis. A muon pair is considered an T(IS') candidate if its invari- 
ant mass is in the signal region defined by 9300 MeV/c^ < M(/x+/x^) < 9600 MeV/c^; 
this selection yields a sample of 2186 events. The number of background events in 
this sample is obtained by fitting the invariant mass distribution to a polynomial plus 
three Gaussians and integrating the function associated with the background in the 
signal region. The resulting number of T(IS') mesons is 1462 ± 55. 

Photon candidates are selected by demanding a transverse energy deposition of 
at least 0.7 GeV in a cell of the central electromagnetic calorimeter and a signal 
in the fiducial volume of the proportional chambers (CES) which are embedded in 
the calorimeter at a depth of six radiation lengths. The fiducial volume requirement 
ensures that the shower is fully contained in a cell. The location of the signal in the 
CES chambers and the event interaction point determine the direction of the photon 
momentum; its magnitude is the energy deposited in the calorimeter. We correct the 
photon energy for the energy lost in the material in front of the calorimeter based 
on a simulation of the detector response to photons. For low energy photons the 



average correction factor varies from 1.03 to 1.14 depending on the polar angle. We 
have verified that the simulation is trustworthy by comparing the simulated electron 
response with the response of electrons from photon conversions found in the data. 

To reduce the combinatorial background resulting from multiple photon candi- 
dates per event we apply the following isolation requirements to the photon: a) 
no charged particle track should point to the photon cell, b) only one CES cluster 
should be associated with the cell, and c) the total electromagnetic energy in the 8 
cells neighboring the photon must be less than 0.5 GeV. The T(IS') is combined with 
all remaining photons within the 90 degree cone around the T(IS'), and the invariant 
mass difference, AM = M(/i"'"/i~7) — M{fi~^fi~), is calculated. The AM distribution 
is shown in Fig. 0. There are two well separated signals; their masses and widths 
are consistent with expectations based on a simulation of the radiative decays of the 
Xfe(lP) and X6(2-P) mesons. The individual angular momentum states of the Xfe's 
(J=l,2) however, cannot be resolved. 

The shape of the background resulting from combinations of the T(IS') with pho- 
tons unassociated with Xb decays is obtained with a Monte Carlo method that uses 
T(IS') candidate events as input. We consider as sources of photons: a) decays of 7r° 
that are not from rj or Kg decays, b) t] decays, c) Kg decays. These sources are simu- 
lated by replacing each charged particle in the event, other than the two muons, with 
a 7r°, 1] or Kg with probabilities proportional to 4:2:1. These proportions follow from 
isospin symmetry and the ratios K^/tc"^ = 0.25 , tj/tt^ = 0.5 |^. Uncertainties in 
these ratios are considered as sources of systematic uncertainty. The response of the 
detector to the photons resulting from the decay of these embedded neutral particles 



is calculated using a Monte Carlo simulation. Applying the Xb reconstruction to these 
events results in a mass distribution that models the shape of the background. This 
model was tested by comparing the Monte Carlo distribution obtained using events in 
the mass sidebands of the T(IS') peak, with the corresponding distribution obtained 
directly from the data where there should be no Xb signal. The two distributions agree 
well as shown in the inset of Fig. ^. The number of Xb signal events is determined 
by fitting the data AM distribution to the sum of the background distribution, with 
an unconstrained normalization, and two Gaussian functions associated with the sig- 
nals. The mass resolution was fixed to the value calculated by the simulation (60 and 
93 MeV/c^). The fit resuhs in 35.3 ±9.0 and 28.5 ± 12.0 signal events for Xfe(l^) and 
X6(2P) respectively. 

The fraction of T(IS') mesons originating from Xb decays is calculated according 
to the equation: 

]\[Xb 

where A^^*" and A^^ are the numbers of reconstructed Xb and T(IS') mesons respec- 
tively, A]- is the probabihty to reconstruct the photon once the T(15') is found, and 
e'^ is the efficiency of the isolation cuts. 

The photon acceptance, A^, is the product of the probability that the photon is 
within the fiducial volume and the reconstruction efficiency of the fiducial photon. 
The geometric acceptance is determined by using a Monte Carlo simulation where 
Xfe's are generated uniformly in pseudorapidity, and with a Pt distribution equal to 
the measured T(15') spectrum ^. The Xb -^ '^{^S) -y decay is generated with a 
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uniform angular distribution in the Xb rest frame. The T(IS') -^ fi^fi^ decay is also 
generated uniformly in the T(IS') rest frame and the trigger simulation is applied 
to the decay muons. Uncertainties associated with the pt spectrum used for the 
production of Xb niesons and with the unknown Xb polarization are considered as 
sources of systematic uncertainty. The photon reconstruction efficiency is obtained 
from the data by applying the photon requirements, except for the isolation cuts, to a 
sample of electrons from photon conversions selected using only tracking information. 
This efficiency is then corrected for the known differences in the detector response 
between photons and electrons. The reconstruction efficiency rises from 17% to 85% 
for photon with Et ranging from 0.7 GeV to 1.4 GeV. For pif > 8.0 GeV/c, the 
photon acceptance is 0.142±0.004(stat) and 0.284±0.006(stat) for Xfe(l^) and Xfe(2P) 
respectively. The large difference is entirely due to the mass difference between the 
parent particles, resulting in different photon energies. 

To study the effect of the isolation cuts we use a Monte Carlo method that uses 
T(IS') candidate events as input. For each event, we generate a vector distributed 
according to the angular distribution of the photon, relative to the T(IS') momentum, 
obtained simulating the decay Xb ~* T(IS') 7. The probability that the isolation 
requirements are satisfied when applied to the calorimeter cell intercepted by the 
vector gives the cut efficiency. Since there are background events in the T(IS') signal 
region, we measure the efficiency in the signal and sideband regions and derive the 
efficiency associated with T(IS') mesons. The resulting efficiency is e'^ = 0.627 ± 
0.013(stat) for Xbi^P) and e^ = 0.651 ± 0M3{stat) for XbC^P); the difference is due 
to the different kinematic of the decays. We assume that this efficiency, calculated 
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from the inclusive sample of T(IS') events, is applicable to the subsample of interest 
where the T(IS') originates from a xt- This assumption is supported by a study using 
samples of J/ip events. We calculate e^ using the inclusive sample of J/ip events, with 
the Monte Carlo method just described, and independently using a pure sample of 
J/ip from Xc decay. The latter is the sample of Xc -^ J /''Pi reconstructed by requiring 
the photon to convert into an electron-positron pair. In this sample we measure the 
efficiency by applying the isolation cuts to the calorimeter cell which would have been 
hit by the photon, had it not converted @]. This measurement yields an efficiency 
of 0.57 ± 0.06(stat); the Monte Carlo calculation is in good agreement, yielding an 
efficiency of 0.56 ± 0.01 (stat). 

The systematic uncertainty on F^j^^^ associated with the Xb production and decay 
model is estimated by varying the shape of the pt spectrum as well as the decay 
angular distribution to account for fully polarized Xb's; the uncertainty is ±13% for 
Xbi^P) and ±9% for X{,(2P). The uncertainty in the determination of A^-^' is ±7% 
for X(,(1P) and ±9% for Xf,(2P). This includes the effect of varying the 7r°, rj and Kg 
composition in our background model from 4:2:1 to all 7r°, and a variation of ±2% 
of the calorimeter energy scale used in the simulation. It also includes the effect of 
varying the resolution of the Gaussians used in the fit by ±6%, the uncertainty on the 
resolution. An uncertainty of ±6% for ^^(IP) and ±3% for Xb(2P) is associated with 
the estimation of the detector response difference between photons and electrons. An 
additional ±4% uncertainty arises from the statistical and systematic uncertainties 
associated with e'^. We combine these uncertainties, assuming they are independent, 
into a total systematic uncertainty of ±16.4% for Xfe(lP) and ±13.7% for X6(2P). 
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The fractions of T(IS') mesons, with p^ > 8.0 GeV/c, that come from Xbi^P) and 
Xb{2P) decays are (27.1 ± 6.9{stat) ± 4.4(sys))% and (10.5 ± 4.4(stot) ± lA{sys))% 
respectively. 

To calculate the fraction of directly produced T(IS') mesons we must estimate 
the fraction of T(15')'s associated with sources other than Xbi^P) and Xfe(2P). We 
calculate the contribution due to T(2S'), T{3S) —>■ T{lS)7ni using a Monte Carlo sim- 
ulation of these decays normalized with the T(2S') and T(3S') cross section measured 
in this experiment 0. We find that the fraction of T(lS')'s, with p^ > 8.0 GeV/c, 
from T{2S) and T(3S) decays, is {I0.7tli)% and (0.8l|]:^)% respectively. An addi- 
tional contribution could be associated to the yet unobserved X6(3-P) mesons. These 
states are predicted to lie below BB threshold and to decay radiatively to T(IS'), 
T(2S') and T{3S). An upper limit on the fraction of T(lS')'s from Xb(3P) decays 
can be calculated with the conservative assumption that all T(3S') mesons in our 
data come from X6(3-P) decays. To estimate the contribution to T(IS'), relative to 
T(3S'), we have used a theoretical calculation of the radiative decay widths of the 
Xb(3P) [Q and the detector simulation to take into account the effect of the trig- 
ger and kinematical cuts. Our estimate is that less than 6% of the T(15')'s, with 
p^ > 8.0 GeV/c, arise from Xfe(3P) decays. We derive the fraction of directly pro- 
duced T(IS') mesons according to the equation F^^-^ — 1 — P^*-^'^-' — Px 7 where 
Py is the fraction of T(15')'s from T(2S') and T(3S'). Systematic uncertainties 
on F^j^l arise from uncertainties on the T(2S') cross section and branching frac- 
tions. Our upper limit on the contribution from X5(3P) decays is also considered a 
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systematic uncertainty, and is added in quadrature to the negative error. We find 
fJS'^^'^ = (50.9 ± 8.2{stat) ± 9.0{sys))% for pj > 8.0 GeV/c. 

In conclusion, we have measured the fraction of T(IS') mesons originating 
from Xb decays and derived the fraction of directly produced T(lS')'s. We find 
that (27.1 ± 6.9{stat) ± 4:A{sys))% of all T{1S) mesons with pj > 8.0 GeV/c 
come from Xbi^P) decays, (10.5 ± 4.4(stat) ± lA{sys))% from Xbi'^P) decays, and 
{50.9 ±8.2{stat) ±9 .0{sys))% are directly produced. A calculation based on the Color 
Singlet Model |T^ predicts a contribution of about 41% and 13% from Xbi^P) and 



Xbi'^P) decays respectively. This measurement will allow the determination of the ma- 
trix elements associated with the production of Xb{^P), Xbi'^P) and T(IS') mesons, 
thus providing information on color-octet contributions in bottomonium production. 
We thank the Fermilab staff and the technical staffs of the participating institu- 
tions for their vital contributions. This work was supported by the U.S. Department 
of Energy and National Science Foundation; the Italian Istituto Nazionale di Fisica 
Nucleare; the Ministry of Education, Science and Culture of Japan; the Natural Sci- 
ences and Engineering Research Council of Canada; the National Science Council of 
the Republic of China; and the A. P. Sloan Foundation. 



14 



[1] F. Abe et al, Phys. Rev. Lett. 79, 572 (1997); ibid. 79, 578 (1997). 

[2] E. Braaten and S. Fleming, Phys. Rev. Lett. 74, 3327 (1995); G.T. Bodwin, E. Braaten 
and G.P. Lepage, Phys. Rev. D51, 1125 (1995). 

[3] F. Abe et al, Phys. Rev. Lett. 75, 4358 (1995). 

[4] P. Cho and A.K. Leibovich, Phys. Rev. D53, 150 (1996); ibid. D53, 6203 (1996). 

[5] F. Abe et al, Nucl. Instr. Meth. Phys. Res., Sect. A 271, 387 (1988). 

[6] F. Abe et al, Phys. Rev. D57, 5382 (1998). 

[7] T. Alexopoulos et al, Phys. Rev. Lett. 64, 991 (1990); M. Banner et al, Z. Phys. 
C27, 329 (1985). 

[8] The smaU conversion probability, and the lower available statistics, prevent us from 
applying this technique directly in the T(IS') sample. 

[9] F. Daghighian and D. Silverman, Phys. Rev. D36, 3401 (1987). 

[10] A.K. Leibovich private communication. The calculation is based on the Color Singlet 
contribution, as described in fH. 



15 



Figure Captions 

FIG. 1. The invariant mass distribution of muon pairs after the selection de- 
scribed in the text. Region S is the T(IS') signal region, region B defines the 
T(IS') sidebands. The solid line is the function used to fit the data, the dotted 
line is the function used to calculate the number of background events in the 
signal region. 

FIG. 2. The mass difference distribution, AM = M(/i+/i^7) — M(yU+/i^), after 
the selection described in the text. The points represent the data. The shaded 
histogram is the background shape predicted by the Monte Carlo calculation. 
The solid line is the fit of the data to two Gaussian functions plus the background 
histogram. The inset shows the comparison between the AM distribution for 
dimuons in the T(IS') sidebands (region B of Fig.l), and the corresponding one 
predicted by the Monte Carlo calculation; the two distributions are normalized 
to equal area and the vertical scale is arbitrary. The size of the bin is the same 
in both figures. 
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